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Deep 610-MHz GMRT observations of the Spitzer extragalactic 
First Look Survey field - 1. Observations, data analysis and source 
catalogue 
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ABSTRACT 

Observations of the Spitzer extragalactic First Look Survey field taken at 610 MHz with the 
Giant Metrewave Radio Telescope are presented. Seven individual pointings were observed, 
covering a total area of ^4 deg 2 with a resolution of 5"8 x A'.'l, PA 60°. The r.m.s. noise at 
the centre of the pointings is between 27 and 30 pi Jy before correction for the GMRT primary 
beam. The techniques used for data reduction and production of a mosaicked image of the 
region are described, and the final mosaic, along with a catalogue of 3944 sources detected 
above ~ 5(7, are presented. The survey complements existing radio and infrared data available 
for this region. 

Key words: surveys - catalogues - radio continuum: galaxies 



1 INTRODUCTION 



The Spitzer First Look Survey was the first major s cientific pro- 
gram carried out by the Spitzer Space Telescope ( Wer ner et al.1 
2004), observing for ~ 1 10 hours of Director's Discretionary Time 
in December 2003. The aim of the extragalactic component of the 
First Look Survey (xFLS) was to study a region with low Galactic 
background to a significantly deeper level than any previous large- 
area infrared survey, in order to accurately characterise the dom- 
inant infrared source populations. Observations were made over 
four square degrees centred on L7 h 18 m 00 s , +59°30'00" (J2000 
coordinates, which are used throughout this pa per) with two in - 
struments - the Infrared Array Camera (IRAC. lFazio et al]|2004l) . 
using the 3.6, 4.5, 5.8 and 8 flm bands, and the Mu ltiband Imag- 
ing Photometer for Spitzer (MIPS. lRieke et alj|2004l) at 24, 70 and 
160 flm. The data have been processed and maps and s ource cata- 
logues are curren tly available (IRA C -lLacv et all 2005; 24 /im - 
iFadda et alj2004 70 and 160 Atm - IFraver et alj|2006t) . 

Complementary observations have been taken at a range of 
wavelengths to fully exploit the new deep infrared data. Deep 
optical surveys of t he region have b een completed in the R- 
band (KPNO 4 m, Fadd a et al l |2004|) . and in the u*- and g- 
bands (CFHT 3.6 m. lShim e t al. 2006). The xFLS region was cov- 
ered by the early dat a release of the Sloan Digital Sky Survey 
( Stoughton et al. 2002). A further redshift survey targeting selected 
24-/j m sources was made w ith the MMT/Hectospec fiber spectro- 
graph (Pap ovich et al.| [2006). and a total of 1587 redshifts are pub- 
licly available. 

There are two existing 1.4-GHz radio surveys of the xFLS re- 
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gion, made with the Very Large Array (VLA) in B-array configura- 
tion, and the Westerbork Synthesis Radio Telescope (WSRT). The 
VLA survey covered ~4 deg 2 wi th a resolution of 5", to a 1 o" depth 
of ~23 /zJy jCondon et al.ll2003h and contains 3565 sources, while 
the WSRT survey jMorgantietai]|2004l) concentrated on ~1 deg 
of the xFLS region to a depth of 8.5 /iJy beam -1 with a resolution 
of 14" x 1 1", PA 0°, and contains 1048 sources. 

There is a tight correlation between the far-infrared and ra- 
dio luminosities of galaxies, whi c h has been k nown for many 
years (see e.g. iHelou et al.l 1 19851: ICo ndon 199 2J). The correla- 
tion applies to both the local JHughes et al.l 120060 and global 
(Murphy etal. 2006) properties of galaxies in the local and dis- 
tant universe Icia rrett 2002; Gruppioni et al. 20031: Ichapman et al.l 
l2005l : lLuo& Wull2005l) . There have been t wo previous studies of 
the IR/radio corre lation in the xFLS region (Appl eton et~ai1l2004l : 
IFraver et al. 2006). Appleton found no evidence for a variation of 
the correlation with redshift, while Frayer, using more sensitive 
infrared data, found a decrease in the infrared/radio flux density 
ratio with z. To study the luminosity variation it is necessary to 
accurately fc-correct the observed radio flux densities to their rest- 
frame values using the spectral index a (where a is here defined so 
that flux density scales with frequency as 5 °= v~ a ). Both previous 
works in this region assumed that all sources have the same radio 
spectral index (taken to be a = 0.7 and 0.8 respectively), but in or- 
der to perform this correction accurately it is necessary to have a 
detection in at least two radio frequencies. 

The large amount of existing data on this field makes it a good 
candidate for a further deep radio survey. We have imaged the 
xFLS region at 610 MHz with the Giant Metrewave Radio Tele- 
scope (GMRT), reaching an r.m.s. noise level of around 30 /iJy 
before primary beam correction. Our survey has a comparable res- 
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Figure 1. The seven pointings observed with the GMRT. Pointing centres 
are spaced by 43 arcmin. The greyscale represents the r.m.s noise level cal- 
culated by Source Extractor (see Section 3), after correcting each pointing 
for the primary beam response of the GMRT and mosaicking the pointings 
together. The 45 jliy contour (equivalent to the 23 £lJy noise level of the 
VLA 1.4 GHz survey, for a source with spectral index of 0.8) has been 
plotted. 
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Figure 2. The uv coverage for the central pointing. Baselines less than 1 kA 
were not used in imaging and have been omitted. 

olution (5." 8 x 4"7 compared with 5") to the VLA 1.4-GHz sur- 
vey, allowing a direct comparison between source properties at the 
two wavelengths to be made. Our observations are deeper than the 
VLA survey over most of the image, for a typical radio source with 
a = 0.8. 

Section 2 describes the observations, and data reduction tech- 
niques. In Section 3 we present the mosaicked image and source 
catalogue. We show a selection of extended objects visible in our 



mosaic, and discuss the artefacts seen near bright sources. In Sec- 
tion 4 we compare our catalogue with the VLA survey, in order to 
test the positional accuracy of our data. Further analyses of the data, 
including detailed comparisions with Spitzer data will be made in 
forthcoming papers. Appendix A presents some technical details 
of corrections made to the observed GMRT data, and Appendix B 
gives details of a correction made to the integrated flux densities of 
sources. 



2 OBSERVATIONS AND DATA REDUCTION 

Observa tions of the xFLS region were made at 610 MHz with the 
GMRT ( Ananthakrishnan 2005), located near Pune, India, over 
four days in March 2004. Seven pointings were observed, in a 
hexagonal grid centred on 17 h 18 m 00 s , +59°30'00" as shown in 
Fig- El The pointings were spaced by 43 arcmin, approximately the 
half-power bandwidth of the GMRT primary beam at 610 MHz, 
which gives nearly uniform noise coverage over most of the xFLS 
region. 

Observations of 3C48 or 3C286 were made at the beginning 
and end of each observing session in order to calibrate the flux 
density scale. The AIPS task SETJY was used to calculate flux 
densities at 610 MH z of 29.4 and 21.1 Jy respectively, using an 
extension of the lBaars~et afl d 19771) scale. Two sidebands were ob- 
served, each of 16 MHz with 128 spectral channels, and a 16.9 s 
integration time was used. Each field was observed for ~200 min 
over four 10 hour observing sessions made between March 23rd 
and 27th. The observations consisted of interleaved 20 min scans 
of each pointing, in order to improve uv coverage. A nearby phase 
calibrator, J1634+627, was observed for three minutes between 
the scans of each pointing to monitor any time-dependant phase 
and amplitude fluctuations of the telescope. The measured phase of 
J 1634+627 varied smoothly between observations, with a typical 
variation between calibrator observations of less than 40° for the 
long baseline antennas and below 10° for the short baseline anten- 
nas. 

Standard AIPS tasks were used to flag bad baselines, anten- 
nas, channels that were suffering from narrow band interference, 
and the first and last 16.9 s integration period of each scan. A band- 
pass correction was applied using the flux calibrators, for each an- 
tenna. A pseudo-continuum channel was then made by combining 
the central ten channels together, and an antenna-based amplitude 
and phase calibration created using the observations of J1634+627. 
This calibration was applied back to the original 128 channel data 
set, which was compressed into 1 1 channels, each containing ten 
of the original spectral channels (so the first and last few channels, 
which tended to be the noisiest, were discarded). The small width of 
these new channels ensures that bandwidth smearing is not a prob- 
lem, and all 1 1 channels could be individually inspected to remove 
additional interference. After the flagging and calibration was com- 
plete the two sidebands were combined into a single data set (see 
Appendix Al) to improve the uv coverage. The coverage for the 
central pointing is shown in Fig. [2] 

The large field of view of the GMRT leads to significant errors 
if the whole field is imaged directly, due to the non-planar nature of 
the sky. To minimise these errors, each pointing was broken down 
into 19 smaller facets which were imaged separately, with a dif- 
ferent assumed phase centre, and then recombined to deal with the 
transformation from planes to the sphere. Images were made with 
an elliptical synthesised beam with size 5l'8 x 4"7, position an- 
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Figure 3. The 610-MHz image of the Spitzer extragalactic First Look Survey field. The greyscale ranges between —0.2 and 1 mJy beam -1 . Artefacts are 
visible around the brightest few sources. 



gle +60°, with a pixel size of l!'5 to ensure that the beam was well 
oversampled. 

The GMRT has a large number of small baselines, due to the 
cluster of 14 antennas in the central 1 km of the array. This domi- 
nates the uv coverage and affects the beam shape, so baselines less 
than 1 kl were omitted in the imaging. The images went through 
three iterations of phase self-calibration at 10, 3 and 1 minute in- 
tervals, and then a final round of self-calibration correcting both 
phase and amplitude errors. The overall amplitude gain was held 
constant in order not to alter the flux density of sources. The initial 
self-calibration step corrected the phase errors by up to 10°, with 
later self-calibration making much smaller changes and having a 
smaller effect on the r.m.s. image noise. 

Two problems were identified and corrected during imaging 
(see Appendix A for further details). 

(i) An error in the time-stamps of the uv data, and hence the uvw 
coordinates, was corrected using a custom-made AIPS task. 



(ii) The GMRT primary beam centre was shifted by ~2.5 ar- 
cmin to compensate for a position-dependant flux density error. 

The theoretical r.m.s. noise of each pointing, before primary 
beam correction, is 



Gsjn{n- 1)N IF AVT 

where T s ~ 92 K is the system temperature, G ~ 0.32 K Jy~' is 
the antenna gain - values taken from the GMRT websit43~ n is the 
number of working antennas, which was typically 28 during our 
observations, A'if = 2 is the number of sidebands, Av = 13.75 MHz 
is the effective bandwidth of each sideband, and T ~ 12000 s is the 
average integration time per pointing. 

The final images have an r.m.s. noise of between 27 and 
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Figure 4. Artefacts visible around a 294 mjy source near 17 h 22 ra 40 s , 
+58°56'00". The greyscale ranges between —0.2 and 1 mJy beam . 

30 /xjy beam -1 before primary beam correction, which is very 
close to the theoretical limit of ~26 /iJy, although dynamic range 
issues limit the quality of the images near bright sources where 
the noise is greater. The seven pointings were corrected for the 
primary beam of the GMRT, taking into account the offset beam 
position as discussed in Appendix A2. The beam correction was 
perfor med using an 8th-orde r polynomial, with coefficients taken 
from lKantharia & Raol OOOlh . The pointings were then mosaicked 
together, weighting the final image by the r.m.s. noises of each in- 
dividual pointing, and the mosaic was cut off at the point where 
the primary beam correction factor dropped to 20% of its central 
value. Figure Q] illustrates the variation in noise across the mosaic. 
The noise level is smooth and around 30 /iJy across the interior of 
the map, and increases towards the edges to about 150 fl]y where 
the primary beam correction was 20%. 



3 RESULTS AND DISCUSSION 

The final mosaic is shown in Fig. [5] There are phase artefacts visi- 
ble around the brightest sources, which have not been successfully 
removed during self-calibration. An enlarged image of one of the 
bright sources is shown in Fig. [4] It is thought that the artefacts are 
due to an elevation-dependant error in the position of the GMRT 
primary beam (see Appendix A2), which will lead to image distor- 
tion near the brightest objects since the observations of each point- 
ing were taken in a series of scans with varying elevations. 

A small portion of the mosaic is shown in Fig. [5] demonstrat- 
ing the quality of the image away from the bright sources, with 
the VLA map of the same area shown in Fig. [6] for comparison. 
The greyscale has been set on Figs [5] and [6] so that an object with 
a spectral index of 0.8 will appear equally bright in both images. 
Most sources are unresolved in the 610-MHz image, although there 
are some objects present with extended structures - we present a 
sample of these in Fig. [7] 

A catalogue of 3944 sour ces was created using Source Ex- 
tractor jBertin & Arnoutsll 19961) with peak brightness greater than 
~ 5<J - see Appendix B for further details. Source Extractor has a 
significant advantage over AIPS tasks such as SAD when creating 
a source list, in that it is capable of calculating the local background 
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Figure 8. Pixel noise distibution, from the Source Extractor r.m.s. map. 
The histogram shows the fraction of pixels at that noise, while the dotted 
line shows the cumulative fraction of pixels with noise below a particular 
value. 
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Figure 9. Source number counts from the 610-MHz Spitzer extragalactic 
First Look Survey field. The solid boxes represent uncorrected counts, while 
the open boxes represent the counts after correction for the survey sensitiv- 
ity level for each solid angle. 

and noise level on the image. Phase errors near the bright sources 
lead to an increase in noise, but by using a box of 16 x 16 pixels to 
estimate the local noise, the number of spurious detections in the 
final catalogue was reduced considerably - this means that the fit is 
less deep near the brightest sources. 

Table Q] presents a sample of 60 entries in the catalogue, 
which is sorted by right ascension. The full table, and radio 
image of the Spitzer extragalactic First Look Survey field will 
be available via http : / /www . mrao . cam . ac . uk/ surveys/ Col- 
umn 1 gives the IAU designation of the source, in the form FLS- 
GMRT Jhhmmss.s+ddmmss, where J represents J2000.0 coordi- 
nates, hhmmss.s represents right ascension in hours, minutes and 
truncated tenths of seconds, and ddmmss represents the declination 
in degrees, arcminutes and truncated arcseconds. Columns 2 and 3 
give the right ascension and declination of the source, calculated 
by first moments of the relevant pixel flux densities to give a cen- 
troid position. Column 4 gives the brightness of the peak pixel in 
each source, in mjy beam -1 , and column 5 gives the local r.m.s. 
noise in [ijy beam -1 . Column 6 gives the integrated flux density in 
mjy, calculated from the mosaic and applying the flux density cor- 
rection factor described in Appendix B. Column 7 gives the error 
in integrated flux density, calculated from the local noise level and 
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Figure 5. A sample section of the GMRT image. The greyscale ranges between —0.2 and 1 mjy beam , and the resolution of the image is 5.8 x 4.7 arcsec - 
PA +60°. Most objects are unresolved, though one clearly extended object is visible near the middle. 
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Figure 6. The existing 1 .4-GHz VLA map for the same region ICondon et all2003l) . The greyscale ranges between —0.1 and 0.5 mjy beam 
the GMRT image for objects with a spectral index of 0.8. The image resolution is 5 arcsec 2 . 
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(i) FLSGMRT J171413.8+594800 and (j) FLSGMRT J172439.4+592044, (k) FLSGMRT J1721 15.4+600432 and 

J171414.6+594750 J172435.7+592053 and J172437.3+592017 J1721 16.7+600406 

Figure 7. A selection of extended objects - contours are plotted at ±100 /ijy beam -1 x 1, 2, 4, 8, . . . Negative contours are represented by dashed lines. The 
resolution of the beam is shown in the bottom left of each image, and the designations o f each source in the GMRT 610-MHz catalogue are given below. All 
sources also appear extended in the VLA 1 .4-GHz catalogue, and images are available in ICondon et all 120031) . 



Deep 610-MHz GMRT observations of the Spitzer extragalactic First Look Survey field - 




17 12 55 50 45 40 35 30 25 20 15 

(1) The extended source has been separated into FLSGMRT J171240.3+595659, J171229.1+595639, J171223.0+595551 and 
J171222.2+595524. The two compact sources are FLSGMRT J171248.1+595528 and J171216.4+595711. 
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(m) FLSGMRT J171 147.7+600836, J171 145.4+600808 and (n) FLSGMRT J171841.3+603630 

FLSGMRT J171 140.8+600731 




(o) The extended object has been separated into FLSGMRT J171647.4+591458, J17 1639.4+ 15 15 1 3, J171631. 1+591544 and 
J171628.4+591527. The double source components are FLSGMRT J171622.4+591536 and J171623.1+591523. 

Figure 7 - continued 
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source size, and taking into account the correction factor. Columns 
8 and 9 give the X, Y pixel coordinates from the mosaic image of 
the source centroid. Column 10 is the Source Extractor deblended 
object flag - for most objects, but 1 when a source has been split 
up into two or more components. For deblended objects it is neces- 
sary to examine the image in order to distinguish between the case 
where two astronomically distinct objects have been split up, and 
when one extended object has been represented by more than one 
entry. There are 211 deblended sources in our catalogue. 

The non-uniform noise characteristics of our survey make it 
important to quantify the area that has been surveyed with each 
noise level. Figure [8] shows the fraction of pixels with a particular 
noise level (taken from the Source Extractor r.m.s. map), and the 
cumulative fraction of pixels at that noise level. The peak bright- 
ness distribution has been plotted in Fig. [9] along with a distribu- 
tion that has been corrected for the varying amount of solid an- 
gle bei ng surveyed to each sensitivity level. Previous studies have 
shown ( Ho pkins et ai] |2002h that Source Extractor has a false de- 
tection rate of below 5%, and detects above 90% of sources with 
peak brightness close to the detection threshold. 
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Figure 10. Source positions for objects detected at 610 MHz and undetected 
at 1.4 GHz. Most sources are located near the centre of one of our pointings, 
where the GMRT sensitivity is greatest. 



4 COMPARISON WITH THE VLA SURVEY 

The VLA survey of the xFLS region l lCondon et alj2003l) has a uni- 
form noise level of 23 /xjy. This is equivalent to a noise of ~ 45 /i Jy 
at 610 MHz, assuming a spectral index of 0.8. Our observations 
have significantly lower noise levels than this in the centre of our 
mosaic, with the noise in the overlap regions between two point- 
ings being approximately 45 /xJy (see Fig.Q3- Figure[8] shows that 
the noise is below 45 flJy for about half the area of our mosaic. 

In order to make a quantitative comparison between the VLA 
xFLS survey and our GMRT survey, we have run Source Extractor 
on our mosaic and the 1.4-GHz image. At the ~5(J level, there 
were 3826 sources detected at 610 MHz, and 3091 detected at 
1.4 GHz, in the region covered by both surveys. The source lists 
were matched using a pairing radius of 6", which is approximately 
the resolution of the two catalogues, and 1580 unique matches were 
found. Fi sure [Tol shows the source distribution of objects found by 
the GMRT but not by the VLA - the majority of the unmatched 
sources are found in regions where our observations are deeper than 
the VLA image. 

It is likely that some of the sources that are not detected at 
610 MHz but are detected at 1.4 GHz are flat spectrum objects, 
with spectral index < 0.5. The VLA noise of 23 /xJy would be 
equivalent to a 610-MHz noise level of ~ 35 flly for a = 0.5, so 
faint flat spectrum objects detectable throughout the VLA image 
would only be detectable in a small region of our mosaic. 

Figure QT] shows the position offsets of the matched GMRT 
sources compared with their VLA counterparts. The offsets have 
an approximately Gaussian distribution, with mean offset in right 
ascension of 0!'4, standard deviation 0"5 and in declination 0"2 
with a standard deviation of 0"6. These offsets have not been ap- 
plied to our catalogue, since it is uncertain as to which survey the 
errors come from. 

The spectral index distribution of objects detected at both 
frequencies is shown in Fig. [T2] The integrated flux densities of 
sources have been used for the calculation, where the flux density 
of VLA sources was corrected using the same method as for the 
GMRT image (details in Appendix B). 
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Figure 11. Source positions in the GMRT catalogue relative to the positions 
found in the VLA image, for unique matches within 6" . 
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Table 1. A sample of 60 entries from the 610-MHz extragalactic First Look Survey catalogue, sorted by right ascension. 
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0.490 
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2371 





FLSGMRT J171751.5+590017 


17: 


17:51.57 
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0.070 


3478 


2724 





FLSGMRT J171753.0+584608 


17: 


17:53.01 


+58:46:08.4 
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FLSGMRT J171755.2+604141 
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0.260 


36 
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0.250 


40 
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0.080 


3439 


1885 






10 



T. Gam et al. 



Figure 12. Radio spectral index between 610 MHz and 1.4 GHz 
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Figure Al. Flux density correction factor for point sources 
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APPENDIX A: CORRECTIONS TO GMRT DATA 
Al Corrections to the uv data 

Two corrections were made to the GMRT uv data before the final 
images were made. 

(i) The uv coverage was increased by combining the two side- 
bands into a single data file. The frequency channels in the lower 
sideband (LSB) have a negative frequency increment, while the up- 
per sideband (USB) has a positive increment - this means that the 
two channels can not be combined directly using standard AIPS 
tasks. We wrote an AIPS task called UVFLP to re-order the LSB 



data, so that it was on the same frequency scale as the USB points, 
and so the sidebands could be combined before imaging. 

(ii) The coordinates of sources originally found in our GMRT 
images were seen to be slightly rotated near the edge of each point- 
ing compared with their VLA positions. This was due to incorrect 
time-stamps being used in the GMRT online software, leading to 
a slight error in the uv data. We wrote a customised AIPS task, 
UVFXT, to increase the time-stamps by 7 s and correctly recom- 
pute the uvw co-ordinates. This was performed before the final im- 
ages were created. 



A2 Primary beam correction 

During the preliminary analysis of the data, we compared the 
properties of several hundred relatively bright sources with peaks 
greater than 1 mjy beam - 1 , visible in the overlapping regions be- 
tween two pointings. This comparison revealed a systematic differ- 
ence between the apparent brightness of sources in adjacent point- 
ings; sources located to the north-west of a pointing were consis- 
tently brighter than the same source when viewed in the south-east 
of an adjacent pointing. We were able to model this as the effective, 
average pointing centre of the telescope being offset by ~ 2.5 ar- 
cmin in a north-west direction, compared with the nominal point- 
ing centre. The amount and direction of the offset was consistent 
between all pairs of pointings. After applying the correction, the 
systematic effects were removed and the r.m.s. flux density errors 
were below 10% near the edge of the primary beam, compared with 
nearer 20% before the correction was applied. 

It is thought that the primary beam offset depends on telescope 
elevation. Our scans of each source were not taken symmetrically 
about source transit, and the detected offset is likely to be an av- 
erage over the individual offsets of each scan. Sources will have a 
slightly different position and flux density in each scan, due to this 
elevation-dependant error. Self-calibration assumes that the source 
remains constant with time, and this is likely to be the reason why 
residual artefacts are seen near the brightest sources. 



APPENDIX B: SOURCE EXTRACTOR 

Source Extractor calculates the flux density of an object by sum- 
ming all pixels greater than some user-defined threshold. For an 
object to be included in our catalogue, we required it to have at 
least five connected pixels with brightness above 2a, and a peak 
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pixel brightness of greater than 5.25(7. Because of the oversam- 
pling of the beam, the peak of a source was taken to be the value 
of the brightest pixel in the island of flux. The peak brightness re- 
quirement for our catalogue was set to the slightly more conserva- 
tive 5.25cr rather than 5o" due to the increased number of spurious 
sources near the edges of the image that were being detected at 
the lower cutoff. 3944 sources were identified and included in the 
catalogue. 

The integrated flux density measured by Source Extractor only 
comes from pixels above 2o". This means that bright sources with 
a large signal-to-noise ratio (SNR) will have almost all their flux 
measured, whereas for faint sources an appreciable fraction of the 
flux is missed. The required correction factor will depend only on 
the peak SNR, for unresolved sources. We modelled a point source 
with a range of SNRs, and varied the exact location of the centre 
of the source across a pixel. Figure lATl shows the correction factor 
that has been applied to an unresolved source in order to obtain the 
true flux density. In our catalogue we give only the corrected flux 
densities - for extended objects that do not have the same shape as 
the beam, it is necessary to examine the image to obtain an accurate 
flux density. 



